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Abstract 

This paper develops the basic analytical theory related to some recently introduced crowd 
dynamics models. Where well posedness was known only locally in time, it is here ex- 
tended to all of M+. The results on the stability with respect to the equations are improved. 
Moreover, here the case of several populations is considered, obtaining the well posedness 
of systems of multi-D non-local conservation laws. The basic analytical tools are provided 
by the classical Kruzkov theory of scalar conservation laws in several space dimensions. 
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1 Introduction 

From a macroscopic point of view, a crowd can be described through its density p, with 
p € L-'^(IR^; M), and assuming that p satisfies a continuity equation of the form 

dtp + diY {pV) = {) . (1.1) 

The vector V is in general a function of the space coordinate x € and of the density p. 
The latter dependence may well be also of functional type since, in general, it is realistic to 
assume that V depends on p through some sort of weighted space average of p. 

A first example of a model of this kind was presented in [61 Section 4]. There, it is assumed 
that pedestrians follow prescribed paths but adjust their speeds to the density they evaluate 
near to their positions. This amounts to postulate a speed law of the form: 

V = v{p*r])v{x) . (1.2) 

The integral curves of the vector field 7 are the trajectories followed by the pedestrians. For 
instance, 'v{x) is the unit tangent at x to the geodesic curve joining x to the destination of 
the pedestrian at x. The convolution p * -q stands for ri{x — ^) p{t, ^) d.^, for a suitable 
non-negative smooth kernel ry. It represents the average density measured, or felt, by the 
pedestrian at time t in position x. 

Below, we extend the results in j6j proving the global in time existence of the solutions 
to (|1.1|) - (|1.2|) . Moreover, we complete the stability estimates with an estimate on the depen- 
dence of the solutions from v and tT, see Theorem 12.21 The resulting model (|l.ip - (|1.2p enjoys 
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further remarkable analytical properties. Indeed, a standard conservation law generates a 
semigroup which is not differentiable with respect to the initial data see [21 Section 1] for 
an explicit example. On the contrary, under suitable conditions, the semigroup generated 
by p.l|) - (|1.2|) turns out to be differentiable with respect to the data. This allows to obtain, 
rigorously, necessary conditions for optimality in various control problems based on (jl.ip - 
()1.2p , see Section [2] below and [H Section 4] . 

Assuming that pedestrians adapt their path to the crowd density they meet lead to the 
model presented in [5j. There, the speed law 

V = v{p){t{x)+l{p)) . (1.3) 

is considered. Again, ~v is the unit vector field describing the preferred paths. But, contrary 
to (|1.2p . here pedestrians may deviate from it, due to the nonlocal term X, which can be 
assumed, for instance, of the form 

i(rt = -.^^k;sL=. (1.4) 

\/l + ||V(p..,)||' 

Again, {p * r]){t,x) is the average density felt by the pedestrian at {t,x), so that ()2.ip states 
that each individual is ready to leave the preferred path in order to avoid regions where the 
crowd density increases. The denominator in ()1.4p is a normalization factor, so that the 
modulus of V is essentially controlled by the function v in (jl.3p . a smooth non increasing 
function that vanishes at the maximal density. We refer to [5j for further justifications of 
the choices leading to ()l.ip - ()1.3p . Below, we show through numerical integrations that the 
formation of lanes, first noted in [5], is present also in the present multi-populations setting. 

For both problems ()l.ip - ()1.2p and (|l.ip - ()1.3p . we then consider the case of several, say 
n, populations. By this, we mean that different groups of pedestrians are considered, distin- 
guished for instance by their destination. This amounts to consider systems of the form 

dtp' + d[vip'V')=0, i = l,...,n (1.5) 

where, in general, depends on the densities of all populations: = V^{pi, . . . ,Pn)- To 
extend all the above well posedness and stability results, we rely here essentially on [3, 
with the improvements in [TJ]. We recall that systems of the form (jl.Sp are considered also 
in [8], where measure theoretic techniques are exploited. The interaction among different 
populations is considered also in [1] through a macroscopic model and in [9] by means of a 
multiscale model. For a general review about crowd dynamics models we refer to 

The next section deals with the n-populations version of ()l.ip - ()1.2p . Then, Section [3] 
is devoted to the analogous extension of (|l.ip - (|1.3p . presenting also a sample numerical 
integration. All proofs are collected in Section [H 

Throughout, we state and prove every result in W^, for a dimension d G N, d > 0, since 
the 2D case contains the same difficulties as the general d-dimensional situation. By we 
denote the interval [0, -|-(X)[. 

2 A Differentiable Model 

The natural generalization of ()1.2p to the case of n populations is 

V' = v\p^ * + . . . + * 7?") v\x) , (2.1) 



2 



so that * 'rf ){x) is an average of the values attamed by p* in i?(x, 1). The map is the 
usual speed law, typically required to be non increasing since at higher densities the mean 
traffic speed is lower. Below, only the regularity of v is used. The vector 'v^{x) is the direction 
of the pedestrian belonging to the i-th population and situated at x G M^. The presence of 
boundaries, obstacles or other geometric constraints can be described through zf*, see for 
instance [3]. 

This section is devoted to the well posedness of ()1.5p - ()2.ip . extending the results in \Ql not 
only for what concerns the number of populations, but also obtaining global in time existence 
and more complete stability estimates. Before stating the main result, we rigorously specify 
what we mean by solution. 

Definition 2.1. Let T > 0. Fix po G L°°(R'^;M"). A weak entropy solution to <fO)-(l2Tl 

on [0, T] is a bounded measurable map p G C'^ ^[0, T]; L^^^^(M'^; M")^ whose i-th component 
for all i = 1, . . . ,n, is a Kruzkov solution to the problem 

dtP' + diy[pW%x))=0 ^^^^^ V%x) = vipHt)*rj' + ... + p-{t)*rj-) TT^x) . 

p\o,x) = pi{x) V ; 

For the definition of Kruzkov solution, see [13] or |1H Paragraph 6.2]. Above, the convolution 
products in the arguments of u * are intended in the sense 

U{t)*rt){x)= f p\t,Ori\x-i)di i = l,...,n. 

The next Theorem summarizes various results in |6], particularized to ()l.ip - ()1.2p . 
Theorem 2.2. Fix d, n € N with d,n > 0. Assume the following conditions: 
(v) G (C^ n W2'°°)(M;M) fori = l,...,n. 
(v) v' G (C2nW2'i)(M'^;Si) /ori = l,...,n. 

(77) rf G (C2nW2'°°)(M'^;[0,l]) and 



= 1 for i = 1, . . . , n. 

Li 



Then, there exists a semigroup 5: M+ x (L^ n L~ n BV)(M'^; M") (L^ n L°° n BV)(M'^; M") 
such that: 

1. For all Po G (L^ n L°° n BV)(IR'^; M"), for all t > 0, the orbit t ^ StPo is the unique 
solution to ()l.ip - (j2.ip in the sense of Definition \2. 1\ with initial datum po- Furthermore, 

the map t ^ Stpo is in C° ('m+; L^(M'^; M")) . 



2. For all po G (L^ n L°° n BV)(M'^; M"), if pi > for i = 1 . . . , n, then {StPo)i > for all 
t > 0. 

3. There exists a constant C such that for all po G (L^ n L°° n BV)(M'^; M"), the corre- 
sponding solution satisfies for all t G M"*" 

TY{p{t)) <{TY{po)+Ct\\poh^)e^' and \\p{t)\\j^^ < \\poho. e^' . 
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4- Fix a positive M. Then there exists functions L, Arj, Ay, Ajf G C'^(M+;]R+) such for 
allpo,i,Po,2 inL^(M°';M") mt/i max |||po_i||j^i, ||po,j||j^^, TV (/9o,j)| < M , for allvi,V2 
satisfying {v), for all iTi, 1^2 satisfying (t?) and for all 771,772 satisfying (rj), the corre- 
sponding solutions pi , p2 satisfy, for all t S M"*" , 

||pi(t) -^2(^)111^1 < {1 + tL{t)) \\po,l - po,2\\j^i 

+ tA 

+tAjrit) (II - ^^2|Il°° + ll^^i - ^^2|lwi.i) • 

5. More regular initial data imply more regular solutions, in the sense that 

Po G (W^'i nL°°)(M'^;R") =^ G M+, p{t) (£ W^'^R'^;W) , 
Po G Wi'°°(R"';R'^) =^ yt£R+, p{t)eW^'°°{ 



Furthermore, there exists a positive constant C such that 

||p(t)||wi,i < (1 + Ct)e^' IIPoIIwi.i and \\p{t)\\^,,^ < (1 + Ct)e^' IIPoIIwi.- • 

6. If V, 'v and rj are of class C^, then 

Po G (W^'i nL°^)(M^;M") =^ G M+, p{t) , 
Po G W2''«(R"';M") =^ VtGM+, p(t) G W2'°°( 
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and for a suitable non-negative constant C , we have the estimates 
\\p{t)\\^,.,<e''\l + Ctf\\ 

Po|lw2.i 1 

\\p{t)\\^,.^<e''\l + Ctf\\ 

Po|lw2.oo • 

7. // V is of class C^ for any initial data po G (W^'"^ n W2'i)(M°'; M"), cjo G (W^'^ n 
L°°)(M'*; M") and for all time t G the semigroup S is strongly L,^ Gateaux differ- 
entiable in the direction Gq- The derivative T>St{po){(^o) of St at po in the direction ctq 
is 

D5t(po)(ao) = Sf°(cTo), 
where is the linear semigroup generated by the Kruzkov solution to 



(2.2) 



dta^ + div (^a'V'ip) + p' DV'{p){a)j =0 (t, x) G / x M'^ 
a'{0,x) = ai{x) xGR"^ 

where p{t) = StPo for all t G and is defined in ()2.ip . 

Note that the hnear problem (j2.2p is the equation that is obtained through a merely formal 
linearization of ()1.5p - ()2.ip . 

The above regularity results allow to state and prove the following necessary condition for 
optimality. 

Proposition 2.3. Let f G Ci'^(M";R+), t/j G L°°(M+ x M'^;M+) and assume that the prob- 
lem lll.l\) - [2A\) satisfies the assumptions at 6. in Theorem \2.2l Denote by S: M^" x (L-*^ fl 
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L°°)(M'^;M") (L^ n L°^)(]R^;M") the semigroup generated by [r^)^ [2l\) . Introduce the 
integral cost functional 

J{po) = I fiStPo) Ht,x)dx . (2.3) 

Then, J is strongly L°° Gateaux differentiable in any direction ao G (W^-^ n L°°)(M'^; M"). 

Moreover, letT.: M+ x (Wi'inL°°)(M'^; M") (Wi'inL°°)(M'^; M") be the linear semigroup 
generated by 112. 2\) . Then, 

DJ{po){cJo)= I f'{StPo)T.P°{ao)il>{t,x)dx. 

If p^ G (L^ nL°°)(M'^;M") solves the problem 

find G (L^ nL~)(M'^;]R") such that J{p^) = min J{p) 

PoG(LinL°°)(Rd;R") 

then, for all Uo G (L^ n L°°)(M'^; R"), 

/ f ' {Stp.) ^'t* ao ip{t, x) dx = 0. (2.4) 

The proof directly follows from [6l Proposition 2.12 and Theorem 4.2] and is hence omitted. 

Remark that (|1.5p - (|2.ip provides an environment where optimal control problems can be 
considered. On the other hand, no uniform upper bound in L°° is available. 

3 Pattern Formation 

In the case of n populations trying to avoid each other, we are led to consider (|1.5p with 

V' = v\p') {ir\x)+r{p\...,p^)) for i = l,...,n, (3.1) 

where are suitable nonlocal functionals. According to (j3.ip . the velocity of 

the i-th population is the product of a scalar crowding factor v^{p^) with a vector lj^[x) + 
X*(/?^, . . . , /j"), which is the sum of a preferred direction 'v"^{x) and a deviation X*(p^, . . . , /o"). 
The scalar v'^{p^) approximately gives the modulus of the speed. A further standard condition 
on typically required in the engineering literature is that t;* be weakly decreasing. However, 
this assumption is here not exploited. The unit vector field tT* can be, for instance, the 
vector tangent at x to the geodesic that the individuals in the i-ih population would follow 
to get to their destination, if unaffected by any other individual. The term Z*(p^, . . . , p") 
describes how the i-th population deviates from its preferred trajectory due to the interaction 
among individuals, both of the same and of different populations. In general, it is a nonlocal 
functional, since its value at any position x depends on the population densities averaged over 
a neighborhood of x. The present setting generalizes the model in [5]. 

Below, we first address the main analytical properties of ()1.5p ~ p.ip . such as the existence 
of solutions, their continuous dependence from the initial data and their stability with respect 
to tT* and X*. Then, a numerical integrations shows further qualitative properties of the 
solutions to (fT3]l -(l3TT). 

Denote by -R > a given maximal density. Throughout, we also denote the flux of the 
i-th population by g*(/0*) = p'^v^{p^), for aU G [0,R]. 

Our starting point is the rigorous definition of solution to (jl.5p - (j3.ip . analogous to Defi- 
nition 12.11 
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Definition 3.1. Fix the initial datum {pi, . . . , p"^) G (L^ n L°°)(M'^; [0, i?]"). A map p G 

C° (^[0,T];L^{R'^;[0,RY'-j^ is a weak entropy solution to (fT3]l - ([311) corresponding to the 

initial condition if, for i = l,...,n, is a Kruzkov solution to the Cauchy 

problem for the scalar conservation law 



dtp^ + div f*(/9*) x) 







where V\t, x) = lj'{x) + V . . . , p„(t)) (x) 



For the definition of Kruzkov solution we refer to |131 Definition 1] . 
Theorem 3.2. Assume the following conditions: 
(v) For i = 1,. . . ,n, i;^ G C2([0, i?]; M+) satisfies v\R) = 0. 
{v) Fori = l,...,n, e (C^ n Wi'~)(M'^; M'^) and div 7* G Wi'i(M'^; R'^^'^) 



(I) There exists a constant C/ > such that the functional : L-^(M'^; [0, iil]'^ 
satisfies, for i = 1, . . . ,n, 



VpGLi(M'^;[0,i?]") 



VP(P) 
Vdiv (^X' (p)) 

Tip) -Tip') 
div 



< C"/ IIpIIli(R'';IR")- 



< Ci\\p-p'\ 



;R")' 



r/ien, t/iere exists a semigroup S : M+ x (LinBV)(R"'; [0,RY'-) (LinBV)(M^; [0,i?]") such 
that 

1. For all po G (L^ n BV)(R^; [0,i?]"), the orbit t ^ Stpo is the unique solution to (jl.ip - 
(|3.ip mt/i initial datum po in the sense of Definition \3.1[ 



2. For allpo G (L^ nBV)(M'^; [0,ii]"), t/ie map t ^ StPo is in C° (^M+;L1(M'^; [0,ii]"; 

5. For all po G (L^ n BV)(M''; [0,ii]"), t/ie following estimate holds: 

TV < TV ipo) e^"' + dWd e"°*|kllLoo([o,,j])(C/ + ||div ir\\j^^)t , (3.2) 

where is defined in (|4.1|) . 

^. Fix M > 0. Let vi,V2 satisfy (v), ^1,^2 satisfy (v) and Ii,l2 satisfy (I). Then, there 
exist b,ce C°(M+;R+) such that for all Po,i,Po,2 G Li(M^; [0,R]'^) with TV (poi) < M 
and for all t G 

Po,l — Po,2||li 

+tcit) (\\qi - g2|lwi.°° + 11^1 - ^2|Il°° + ||div(iri - tr2)||Li) • 



and the functions b,c depend ond, Cj, R, \\qi\\^ffi,oo, || i|lwi.°°' ll^iv i^i 11-^^^1,1 , ||Vr/||Lo 
and on M . 
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Detailed expressions for the functions b and c are available, together with the proof, in Para- 
graph [Ol 

This theorem also provides a kind of maximum property since each density remains 
bounded by R. However, it does not guarantee that the sum pi{t) + . . . + Pn{t) remains 
bounded by R. 

Preliminary to any use of ()1.5p - ()3.ip is the choice of a specific X. First, we consider the 
following lemma that eases the construction of operators that satisfy (I). 

Lemma 3.3. Let N:W^ he defined by 

11 

N{u) = 



1 , II l|2 

1 + \\u\\ 



If I satisfies (I), then so does MoT. 

The proof consists in long but elementary computations and is hence omitted. By 
Lemma [3. 3 1 it is immediate to check that ()1.4p is satisfied by the operator X defined in ()1.4p in 
the case of one population, as soon as r/ G C^(M'^; M+) and = 1, considered in [6j. When 

more populations are present, it is natural to consider different kinds of interactions. For in- 
stance, the population p^ might deviate from its preferred path 'v^ due to the population p^ 
pushing in the direction tT^, thus leading to consider the operator 

^ (p) = I ^ ^ ' , (3.3) 



1 + 



'^p'^v{p^)lij'^^ * rj 
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Under assumption (v) on v and with rj G C^(]R'^; R"*") such that = 1, also as defined 

in ()3.3p satisfies (I). 

When these or other operators are to be considered together, then the following lemma 
makes the verification of (I) immediate. 

Lemma 3.4. Let X and X satisfy (I). Fix any two functions a,l3 £ (C^ n W^'°°)(M'^; M). 
Then, also aX + f3I satisfies (I). 

The proof is elementary and hence omitted. 

Differently from the model in Section [21 the semigroup generated by ()1.5p - ()3.ip is not 
proved to be differentiable with respect to the initial data. On other hand, it develops 
solutions with an apparently rich structure. Moreover, as stated in Theorem 13.21 an upper 
bound in L°° is available. 



3.1 Numerical Integration 

The study of the qualitative properties of the solutions to ()1.5p ~ ()3.ip is in general not amenable 
to purely analytical tools. Numerical integrations, besides being of interest from the point 
of view of the applications, show interesting pattern formations. In the case (|l.ip - (|1.3p of a 
single population, the formation of queues was already noted in [5j. 

The algorithm used is the Lax-Friedrichs method with dimensional splitting. A uniform 
grid {xi,yj) ioi i = 1, ... , and j = 1, . . . , is introduced and the density p is approximated 
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through the values pij on this grid. At every time step, the convolution in T{p) is computed 
through products of the type Ai^ph^B^j^ where the matrices A and B depend only on r/. 

In both the integrations, we use the same geometry and the same numerical parameters. 
More precisely, the space available to the pedestrians is the rectangle M x [—3,3], while the 
numerical domain is [—8,8] x [—4,4]. Pedestrian may exit along the segments {—8} x [—3,3] 
and {8} x [-3, 3]. The mesh size is Ax = Ay = 0.025. 

The preferred path of each pedestrian is the sum g + 5. The vector g is tangent to 
the geodesic towards the pedestrian's target or when the pedestrians would not move. 
The vector 6 describes the discomfort of pedestrian when walking too near to a wall. It is 
perpendicular to the walls and pointing towards the interior of the room. Analytically, it also 
ensures the invariance of the room, see [5j for more details. Its maximum modulus is (5max 
along the walls at |rE2| = A and decreases linearly towards the room interior, vanishing at a 
distance 5r from the walls. In the present integration we set 

<5max = 0.8 , 5r = 0.75 . (3.4) 



3.2 Two Groups of People Crossing 



A classical situation considered in the engineering literature is that of two groups of people 
moving in opposite directions and crossing each other. Typically, lanes, also called paths 
or trails in the engineering literature, are formed, see for instance [12^ ITO] or [1] for a one 
dimensional description. They consist of people going in the same direction. The model (jl.ip - 
()2.ip captures this phenomenon. Indeed, we consider ()1.5p - ()2.ip in the following setting: 







dtp^ + div 

< 




dtp^ + div 


P'v{p^){^ 







yi+||V(pi*»7)| 
yi+||V(pi*r,)| 



£12 



^22 



V(p2*2) 



yi+||V(p2*^)| 
yi+||V(p2*^)| 



(3.5) 



where the p^ population moves to the right and the p^ population move to the left. Above, 
we assume that each pedestrians wants to avoid entering regions with increasing densities and 
that the repulsion towards pedestrians of the other population is bigger than that due to the 
proper population. For the sake of simplicity, we use the same kernel r/ and the same speed 
V in both equations. We consider the specific situation defined by 







+ 5, r]{xi,X2) 
+ (5, v[p) 



[l-(2xO^]'[l-(2x2)^]^X[_„.,_„.,j,(xi,X2), 
4(1-P), 



^2l3, 



(3.6) 



£21 



0.3, £12 
0.7, £22 



0.7, 
0.3. 



As initial datum we choose 

pl{xi,X2) 
pl{xi,X2) 



0.9 ■ Y, , , Sxi,X2) 

'^[-6.4,-3.2] X [-2.4,2.4]^ ' ^' 

0.7 ■ Y, , , Xxi,X2) ■ 

'^[3.2,6.4] X [-2.4,2.4]^ ^' 



(3.7) 
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The resulting numerical integration in Figure [T] shows several lanes forming. First, when most 
of the populations are still separated, the lanes in the two populations are rather symmetric, 
in spite of the different initial densities Then, when the interaction between and gets 

^ Time= 0.000 ^ Time = 1.466 , Time = 1.971 , Time = 2.981 

3 3 

-3 -3 

" -8 -6 -4 -2 2 4 6 8 ' -8 -6 -4 -2 2 4 6 8 "-8 -6 -4 -2 2 4 6 8 ~-8 -6 -4 -2 2 4 S 8 

Time = 0.000 ^ Time = 1.466 , Time = 1.971 , Time = 2.981 

3 3 

2 - I I - 2 

- - 

-2 - - -2 
■3 -3 



Figure 1: Numerical integration of (j3.5p with data (j3.7p and parameters as in (j3.4p - ()3.6p . 
Above, the population moving to the right and, below, p^ moving to the left. Note the 
lanes that are formed. First, due to the self-interaction within each populations and then due 
to the crossing between the two populations. The latter lanes of different populations do not 
superimpose. 

relevant, lanes sharply bend their trajectory so that they are not superimposed. In this way, 
pedestrians lower how much they block each other. After the interaction, the density is so 
low that no more patterns arise. 





3.3 Part of an Audience Leaving a Room 

A sample situations developing various interesting features is the following. Two populations 
are initially uniformly distributed in the same region. At time t = 0, the first populations 
starts moving towards an exit, on the right in Figure while the second moves only to let 
the first one pass. We describe this situation through ()1.5p - ()3.ip . The preferred path of p^ is 
tangent to IT^ = g + (5, where g is the geodesic vector field directed toward the right exit and 
5 is the discomfort, as above. The p^ population then deviates from this trajectory only to 
avoid the other population. On the contrary, the preferred path of p^ is tangent to tT^ = 5, 
since this population moves only to avoid the walls and the p^ population. All this leads to 
the system 



dtp^ + div p^ v{p^) 'v^{x) 



dtp^ + div p^ v{p^) 'v'^{x) 



We consider the specific situation defined by 



'1+ V(p2*,j) 



0. 



0. 



(3.8) 



1+ V(pi*-»y) 



+ ,5, 7?(X1,X2) = [l-(2xi)2]=^[l-(2x2)2]'x._o.5o,5p(xi,X2), 



213, 



(3.9) 



v{p) 



4(1-/,) 



0.3. 
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with initial data 

= 5 y 

^° ' "^[-6.4, -3.2] X [-2.4, 2.4] ' 

The resulting numerical integration is in Figure El Note that the initial distributions of 



1,2, 



(3.10) 



Time = 0.000 



Time= 0.203 



Time= 1.213 



Time = 2.476 




Time = 0.000 



Time = 0.203 



Time= 1.213 



Time = 2.476 



-8 -6 -4 -2 2 4 E 8 -8 -6 -4 -2 2 4 6 -8 -6 -4 -2 2 4 6 -8 6 -4 -2 6 2 4 8 8 

Figure 2: Numerical integration of (j3.8p with data and parameters as in p.4p - ()3.9p . Above, 
the population and, below, f? . Note first the formation of small clusters separating the 
two populations, then lanes and, finally, a sort of fingering. 

both populations are uniform. Very quickly, patterns start to form. First, clustering: the 
different populations separate forming separated peaks of high density. Then, population 
starts moving significantly to the right along lanes that follow paths in regions of low 
density. Finally, the two populations are almost separated, with a sort oi fingering remaining 
in the region when they are superimposed. Concerning the initial clustering, we underline 
the analogy with [9], where a population is described through a measure possibly containing 
also an atomic part. Concerning the latter fingering, we borrowed here this term from the 
entirely different case of, for instance, thin films, see [15j). 



4 Technical Details 



Below, for any /) G 



pIIl1(K'*;M") stands for J27=i 



. Similarly, TV (/)) 



^^]^ TV (/>*). More generally, if not otherwise stated, the norm of a vector v in M" is the 



1-norm, i.e. \\v\ 



Era 
i=l 



The following constant enters several estimates below: 

»7r/2 



Wd = / (cose)'^d0 . 
Jo 



(4.1) 



The analytical tools below are based on the classical Kruzkov work [13], see also |1H 
Chapter VI]. The stability estimates are taken from [71 114]. 

First, we study the following Cauchy problem for a scalar non-linear conservation law: 



dtp + d[v{q{p)V{t,x)) = 
P(0) = Po . 



(4.2) 
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V{t) G W^'°° 



for all t € 



Lemma 4.1. Assume that 

q € C2(M+;M+) satisfies q{0) = 0, 
V € C^{R+ xR<^;R'^) satisfies 
Po G (Li nL°°)(M'^;M+). 

Then, there exists a unique weak entropy solution p E ^M+, L"'^(M'^; M^)"! to 
If furthermore po € BV( 
= ||p(0||L°o([o,T]xRd); 



(4.3) 



), then p{t) € BV(M'^;R) for all t G M+ and, denoting 



TV(p(t)) <TV(po)e'^°* + Qe'^°*||g||L^(ro^ ,) T / || V div ^(t, dx dr , (4.4) 

and Cd = dWd, with Wd defined in ()4.4p . 



where Ko = {2d+l)\\q'\\^^^^Q^j^^^^\\VV 



'II' n L,— (lUj-tlTl) " i-L°° ([0,r] Xii-i ; - _ „ 

Let now vi,V2, Vi,V2 and Po,i,Po,2 satisfy (|4.3p . Call pi,p2 the solutions to 



dtpi + div {pi vi{pi) Vi{t,x)) = 
Pi(0) = po,l 



dtP2 + dw{p2V2{p2)V2{t,x)) =0 
P2(0) = po,2 



Then, renaming Rt = max{||pi||Loo([o,T]xRd)' IIP2|Il=°{[o,t]xR'«)}' 

||/5l(i) - P2(i)||Li ^ IIPo.I - Po,2||li 



+te'^°* (lk2||Loo([o,R^])ll^l - ^2|Il- + U - l2Loo(^[0,R^]^\\Vlho 

TV(po,i) + Cd l|gi|lL°°([o,i?T]) ||Vdiv Vi(t,x)|| dxdr 



+ lkillL-([o,RT]) J 

+ \\Q1 - 92|1l°o([0,_Rt]) 

Proof. We consider the following steps separately. 



div(Vl(T,x)-y2(T,x)) 

/ / |div V2(r, x)| dxdr 
Jo Jr'^ 



dx dr 



Existence. Let f{t,x,p) = q{p)V{t,x). Then, we have / € C°(M+ x M'' x M;M'^) and 
fit, • , • ) ^ C^(R'^ X ]R;IR'^) for any t G M"^. Through elementary computations the following 
implications can be checked: 

dpf{t,x,p) = q'{p)V{t,x) =^ dpf is bounded on [0, T] x M'^ x [- A, yl] 
div flt,x,p) = q{p)dwV{t,x) div/ G L°°([0, T] x M"^ x A], M) 

dpdiv flt,x,p) = q'{p) diYV{t,x) dpdivf G L'^{[0,T] xR'^ x[-A,A]) 

for all A G M'*'. Hence, (ISj Theorem 4] can be applied and the existence of solutions follows. 



Maximum principle. Since q{0) = 0, p = is solution to (j4.'2|) . The maximum principle 
of Kruzkov [13, Theorem 3] then ensures that Po ^ implies p{t) > for all t >0. 
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BV bound. The BV bound follows from [14:\ Theorem 2.2], which can be applied since 
\7dpf = q'{p) VV{t, x), so that Vdpf G L°°([0, T] x M'' x [-A, A]). Moreover, for any A>0, 



T 



|Vdiv/(t,x, •)||l-'([-a,a])^^'^*= lk(^)|lL-'([-^,A]) y^J|Vdiv ir(t,x)|| dxdt<+oo. 
Let Ko = (2d + l)||g'||Loc([o,i?y])l|V^llL°°([o,r]xMd)- Then, for any i > 0, we have: 

TV(p(t)) <TV(p„)e'^°* + Cde''°*||g||Lo.(ro^ 1) t I ||Vdiv F(t,x)|| dxdr , 
with Cd = dWd, and Wd as defined in (|4.ip . 

Stability estimate. The stability estimate follows from [14i Theorem 2.6]. We have, for 
any ^ > 0, 



||div(/i - /2)(t,x, •)||l-([-a,a]) 



dx dt 



= lki(Pi)||L°°r[-AAn / / |divVi(T,x) - divy2(T,x)| dxdt 

^' ' " Jo JIRd 

+ hi - Q2\\i,oo(!AA\) / \divV2{T,x)\dxdt<+oo. 
' " Jo JlRd 

Then, with Rt = max{||pi||L.c([o,T]xRd)' ll/'2|lLoo([o,T]xRd)}, we get 

\\plit) - P2{t)\\j^i < \\po,l - Po,2\\j^i 

+te-"' (lk2Loo([o,^^])l|V^l -^2|Il- + Iki -'?2||l°«([o,h^]) 



TV(po,i) + Q||gi||Lo 



+ lk2|lL-([0,RTl) / / div(¥L(T,x) -y2(T,x)) 
Jo JM'* 

+ -92|lL°o(ro,ijyl) / / |divyi(T,x)| dxdr 

JO Jw^ 



V div (r, x) 1 1 dx dr 
dx dr 



completing the proof. 



□ 



Corollary 4.2. In the same setting as Lemma \4-1\ if q{p) = Qiip) = Q2{p) = P then q'{p) = 1, 
q"{p) = and the stability estimate reduces to 

\\pl{t) - p2{t)Li < \\po,l-po,2 



iLl 

+te''°*||yi - V, 



+ \\P2\\ 



1 ~ 

TV(p,,i) + C,||pi|| 

L°°{[0,T]xIR' 



L-([o,T]xR'*) / ||Vdiv yi(r,x)|| dxdr 




div (yi(T,x)-V2(T,x)) 



dxdr . 
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4.1 Proofs related to Section [2] 



Proof of Theorem [221 1- Let po G (L^ n L°° n BV)(M^; M") and let A^i = ||po|Ili- Let us 
introduce a given time T > and the sphere 



1 tmd. wnw 



for ah t G [0,T] , ||p(t)|L < iVi L 



equipped with the distance induced by the norm || • |Il°o([o t] li(ir''-ir"))- Consider first a fixed 
i G {l,...,n}. Choose any ri,r2 G Xn^ and denote p\, p\ G L°° (^[0, T]; L1(M'^; R)) the 
solutions of the Cauchy problems (j4.2|) with respectively, for k G {1, 2}, 

qk{p)=P. Vk{t,x)=v\rl*ri^ + ...+rl*'n^). (4.6) 

Thanks to the properties of f * G and r/ the hypotheses on Vk in Lemma 14.11 are satisfied. 
Hence, we get existence and uniqueness of a weak entropy solution, Besides, the bound on 
the solution = ||Po|Ili ^^^Y * ^ ensures that the solution of (|4.2|) - (|4.6p belongs 

to AfjVi • 

Noting furthermore that 

V' G L°°([0,r] X M"';M'^) , 
VV' G L°°([0,r] xM"';M'^^'^)nL'^ ('[0,r];Li(M'*;M'^^^; 



we can apply ^ Corollary 5.2 and Lemma 5.3] or more directly |14^ Theorem 2.2 and Theo- 
rem 2.6]. 

Remark that A'^i > sup^gjo j-] ||^(i)||Li) so that we have the estimate 



div V\t) 



< 



< 





+ 














Loo 





w 



,^(A^i!|V??||loo + 1). 



Thus, with 

Ki = II 

we obtain by [6l Corollary 5.2] 

Kit 



(4.7) 



P\t) 



< 



Loo 



Po 



Loo 



Summing over i = 1, . . . , n, we obtain 



< IIpoI 



We want now to apply Lemma [4. II or Theorem 2.2]. Note that, with the notations of [14j . 
we have 



(2d + l) 



Loo 



< {2d+l)Ki, 
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V div V 



< 



Loo 



El 



LI 



+ 

+2 



y\\A\ 



Loo 



+ 








LI 




Loo 





LI 



< 



IW2.: 



Wl|Vry||^,oo+2iVi||V77||wi.oo + 



Denoting 



we obtain 



K2 = Ibllwi.- Il^llwi.i ^1 !|Vr/||Loo +2A^i !|Vr/||wi,oo + 1 , 



(4.8) 



TV j < TV (pj,)e(2^+i)^i* + te^^^dWdK2 

Let us now apply Corollary 14.21 or [Hj Theorem 2.6]. We obtain 



P1-P2 
where 



< Kte^*||ri -ral 



K 



L°°([0,T];Ll) 



TV (pj,) + (t dVFrfKs + 2 + A^i||Vr?||Lo 











LOO 



M II^I1l°°' ll'^llw^.i) ll^llw2.°o; 



Note that Ki, K2 and K are constants depending on ||v||-y^2,o, 
Ni and d and not on the initial condition po- We do not need here precision on the con 
stant so, with a C large enough, not depending on Pq and X", we have, denoting F(t^ — 
C t e*^* (TV (po) + C{t + 1) ||/Oo|Il°°) ™d summing over i = 1, . . . , n, 



IIpWIIloo < \\po\\-Lc.e^\ 
TV {pit)) < TV(/5o)e^* + C7te^*||/5o||Loo, 



L°°([0,T];Li) • 



Choose T so that F{T) = 1/2. Then, the map Q: -Yat^ — )• ^at^, defined by Q(r) = p, is a 
contraction. Banach Fixed Point Theorem ensures local in time existence and uniqueness of 
weak entropy solutions to (|1.1|) - (|2.1|) . 

In order to get global in time existence, we iterate the previous procedure. Starting from 
time Tn, we obtain 

Upi - P2)(i)||Li < lln - ^2|lLoo([o,t];Li)C(i - r„)e^(*-^") 

X (tV (po)e^^" + CT„e^^"||po||Loo + C(l + t - T„) e^^"\\poho 

Iteratively, we choose Tn+i > Tn such that 

C [Tn+i - Tn) e^^"+i (tV [po) + C (1 + r„+i)) IIPoIIloo) = \ ■ 
Since the sequence Tn grows to +00, we have global in time existence, proving point 1. 
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2. The positivity of the solution directly follows from the fact that />* = is a solution 
and from Kruzkov Maximum Principle [T^l Theorem 3]. 

3. To prove these estimates, we apply [6l Corollary 5.2] and |14^ Theorem 2.2 and Theo- 
rem 2.6] to each equation of the system, which is possible since the coupling is only present 

in the nonlocal term. Let i £ {1, . . . ,n}. Denote Ki as in ()4.7p . A''i = ||/Oo||li = Yl^=i 

and = v^{p^ * r]^ + . . . + * rf") we have 



|p(i)||TOO < \\P0\ 



Loot 



To prove the estimate on the total variation, apply [141 Theorem 2.2]. Denoting K2 as 
in (|4.8p . we obtain 

TV (p\t)] < e^^'i+i)^!* TV (pi) + te^^'^+^^^^'dWdK2 



Summing over i = 1, . . . , n we obtain the desired bound on TV {p{t)) = Yl^=i TV yp^{t)j . 

4. To prove the stability of the solution with respect to initial conditions and param- 
eters, we apply [HI Theorem 2.10]. Denote by p]^, for k = 1,2, the solutions to the Cauchy 
problems 



where Vi{t, x) = v^pl * ril + . . . + pi * tj^) TT^x) . 

for k = 1,2, so that the necessary hypotheses 



dtpl+di^[^3lVi{t,x))=Q, 

Note that divl4* € ([0,r];Li 
hold. Moreover, 

diY{Vl-Vi){t) 



Li 



|Po,i||liI|V7?i||loo + VI-V2 

+ ||po,2||liII^?i - 

|Po,i||liII^^i|Il°° + di"^ 

+ ( ||(Pl - P2)(i)||Lil|VmllL- + ||/5o,2||liI|V7?1 - V7?2|Il° 



+ {\\ipl - P2)(i)||Lill'/l|lL° 



div v'l 



iv\y 



Li 



Li 



+ div {v'l - v'l) 



V 1 — V ' 



Li 



|Po,2||Lil|Vr/2|lL°o , 



and 



< \v\-vl 

+ {Wipi - p2m\ 

Let us introduce 







+ 






— *2 — *2 












LOO 


LOO 


LOO 



+ ||Po,2||Lillm - ??2|Il°° M) 



Loo 



a 



hlilL-ll'?lllL°°||Po,l|lLill^^l|lL-ll'^lllLi + ||4||l°oII^^iIIl°°II'^iIIli 
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+lhi||Locll'7i|lL-l|div iTiIIli , 
/3 = II<IIt.oc||Po,2|| Li|ko.l|lLill^^l|lL°°ll^l|lLi + ||'"l||Li|ko.2||Lill^^'^^lllLi 
+ lk2||Loc||Po,2||Li||^^l|lLi , 

7 = ||div iTiIIli + ||po,i||Li||V?7i||Loc||tri||Li , 

^ = ||^2||lo°I|/'o,2||li||V772||l<x. + |b2||L<x. , 



and 



= Ik'lllT.ooll^lIlL- ' ^' = |hl||L<-lko.2||Li||'yi||Loo , 7' = II^^iIIl- ' = 11^2 Hl- , 



these coefficients being chosen so that 



+5\\v'i - Ij 



2|lwi.l ' 



in-yi){r) < a'\\{p,-p2){T)\\^,+P'\\vi-V2Lo.+j'\\vi-V2\\j^. 



+5'\\'vi - 'V2\\l,oo ■ 



Hence, we get 

pi,i - pi,2 

+ max 



< 



Li 



Li 



TV{po,i) + tdWdK2\\po,i\ 



iVl-Vi){r) 



dr 



{lklW|lLo-||P2(t)|Loo}^ \\dw{V^~Vi){ 



Li 



dr 



< 



Po,l - Po,2 



ip\-pl)ir) ^^dT + tA,{t)\\vi-V2\\^i 



+t Ar^{t)\\r]i -?72||wi.°° +tAt{t) {\\v'i - l^2||wi.i + ll^l ~ ^2||l°c) > 
where Ap, A^, A^i, Ajf are smooth, positive, and increasing, functions of t depending on 
d, ||f||w2.oo, ||i^ilL=°> lhllwi.°°' TV(po,i), max|||po,i||L<x., ||Po,2||l<x.}, ||Po,i||li and 

||Po,2||li- More precisely, denoting 

f(t) = e(2<^+i)^i 
R 
K 



TV(po,i) + idW^d^2||Po,i||L»; > 

max(||Po,i||j^^, ||Po,2||loo) , 



max 
fce{i,2} 



{||i^fc|lwi.oo||l^A;|lwi.-(||Po,fe||Lill^%llL-= + 1)} ' 



we have 

Ap{t) = a'f{t) + aRe^^ , = P' f{t) + pUe^^ , 

Av{t) = 77(0 + 7-^6^^*, A^{t) = 6'f{t) + 6Re''K 
We conclude applying the Gronwall Lemma, obtaining 



(Pl - P2)(*)||li < (||po,l - Po,2||li + tA^{t)\\r]i - ?72||wi.- + tAit)\\vi - V2\\y^i,, 
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+tAjr{t) {\\v'l - 72||loo + II t^i - 'V2\\^fl;l,l 



X ( l + t{f{t)a' + Re^^a) e 



and the stability estimate at 4. follows. 

5. These properties follow through standard computations from the representation for- 
mula given in [6l Lemma 5.1]. 

6. As above, we deduce this further regularity property and the W^'^ estimate from the 
representation formula in [6l Lemma 5.1]. 

7. Let p G C° (^M+; L^(]R'^;M")^ be a solution to the initial problem (fT3]) - (f2T^ such that 
for all t > 0, p{t) G W2'°° n W2'i(M'^; M"). Consider now the linearized equations 



0. 



0. 



where p = {p^ , . . . , p"), r] = {-q^, . . . , r/"), a = {a^, . . . u") and p-k-q = p^ *rj^ + . . . + p"^ * 77". To 
prove the existence and uniqueness of weak entropy solutions to this linearized problem, we 
use the technique that proved to be effective for the initial value problem: let Mi = ||<7o||li, 
let si, S2 € are fixed functions; we fix the nonlocal term and study the Cauchy problems 



dtal + div ( (7* v''-{p-k ri)v^{x] 



div 



for 



o,k 



Or, 



k 

i 



1,2 

l,...,n 



We study the map T: Xmi — > ^^Mi defined by T(s) = a. This application is well defined 
thanks to Kruzkov Theorem [13] and thanks to [6l Lemma 5.1]. Indeed, denoting 



f{t,x,u) 
g{t,x,u) 



uv^{p-krj) iT^^x) 
uv^[p-krj) 'v'^{x) 



F{t, X, u) 
G{t, X, u) 



div ( (f *)'(p -kf]) si-k r]lf^{x) 1 , 
div (p^ {v'^y{p'kr]) 82* r]lj'^{x)] , 



we have duf,dug £ L°°{[0,T]xR'^ x[-U,U]), F-div /, G-div 5 G L°°([0, T] x M'^ x [-[/, [/]), 
duiF - dwf),du{G - div 5) G L°°([0,r] x M'^ x [-U,U]). Hence Kruzkov hypotheses are 
satisfied. To apply [141, Theorem 2.6], we have now to check that 

Vduf €L^{[0,T] xR<^x [-U,U]), 



/o /M^I|V(i^-div/)| 
/o'/M^II^-G-div(/-5)| 



L°°{[-C/,(7]) 
L°°{[-C/,(7]) 



dx dt < +00 , 
dx dt < +00 . 



We have 



\VFit,-,u)L, < \\pit)\ 



W2 



Iw 



1 .00 

II^iIIli II'?IIw2>oo I 



9 + 6||p(t)||Li||V7?||Loo + ||p(t)||Li||Vr?||^. 



+\\pit)Li 
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||Vdiv/(t, •,m)||j^i < \u\K2, 
\\{F-G)it,;u)\\^, < \\p{t)\\ 

X||(S1 - S2){t)\\^, , 

||div(/ - fi-)!*, •,^i)||Li = 0' 
\\dM -9Kt,;u)\\^^ = 0. 

where K2 is defined as in ()4.8p . Using the estimate on ||/o(t)||-^i,oo obtained in 4. and denoting 



K4 — l|po|lwi>°° ll^'llwi'°° ll^llwi.°° II llwi'i (2 + IIPoIIli II^^IIl°°) 



(4.9) 



we obtain 



< 



div (^p' (v^Yip -k r]) (si - S2) * ?? tT^) 



dt 



< ti^4(l + Ct)e'^*||si-S2||L 



= ([0,t],Li) 



For T small enough, we obtain that T is a contraction, proving the local in time existence 
and uniqueness of solutions by Banach Fixed Point Theorem. We then extend to +00 the 
time of existence by iteration of the process. 

Denote now by p, respectively ph, the solution to the initial problem ()1.5p - ()2.ip with 
initial conditions po, respectively po + hao- Moreover, call a the solution to the linearized 
equation (|2.2p with initial condition ao and define Zh = p + ha. If o" is smooth enough we 
can write for the equation 



dtzl + div [ zllv'-ip-kr]) + h{v^) {p * ??) a ★ 77 ) iT* (x) = /i div ( a* (v*) (/> * 77) cr * ri'v\x) 



We want now to estimate \\{ph — Zh){t)^j^^/h. In order to do that, we use [HI Theorem 2.6]. 
As it is similar to the estimate in the proof of [6l, Theorem 2.10], we omit it. □ 

4.2 Proofs related to Section [3] 

Lemma 4.3. In the same setting as Lemma \4-l\ if furthermore for a given R > 0, we have 
q(R) = 0, then 

Po e[0,R] =^ Vt > , p{t) G [0, R] . 

Proof. Since q{0) = q{R) = 0, p = and p = R are solutions to (|4.2p . The maximum 
principle of Kruzkov |13j then ensures that < p < R implies < p{t) < R for all t > 0. □ 

Proof of Theorem 13.21 Consider the following steps separately. 



Existence. Fix an arbitrary positive time T, whose precise value will be chosen later. Let 
(r\...,r'^),(s\...,s") G C° ([0, T]; L1(M'^; [0, i?]")) . For any i G {l,...,n}, we define 



V\t, x) = ir\x) + T r\t), r'^it) (x) , W\t, x) = iT'ix) + V s\t), (x) 
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As V satisfies (v) and I satisfies (I), for any r, s G C° (^[0,T],L1(]R'^, [0,i?]")j, we have for 
any t > 0, V{t),W{t) £ (C^ n Wi'~)(M'^, M"^) and div V (t) , div W (t) G Wi'i(M'^,M). Then, 



according to Lemma [4 .11 there exists p,(J £ C ^[0, T]; L^(M'^; ) ) , that are the weak entropy 
solutions to the systems of decoupled equations 



dta^ +diY {a^v^{a^)W^{t,x)) =0, 
dta"^ + div (fT"t;"(fT") x)) = ^ 



dtp^ + div [p^v^{p^) V^{t, x)) = , 
+ div (p"v"(/9") x)) = , 



with initial condition p(0) = a(0) = /Oo G (L^ n L°°)(M'^; [0, Rf). 

Note furthermore that q^{R) = 0; thus, thanks to Lemma [4.31 that if /? is a solution of the 
above equation, then po G [0,R] implies for all t>0, p{t) G [0,i?]. Hence, we have invariance 
of the interval [0, i?]. 

Consider the map 



T: < 



CO (^[0,r];L^(M'^;[0,i?]"; 
r 



CO (^[0,r];Li(M'^;[0,i?]") 
^ P 



Use Lemma l4.ll and define = maxj < {2d + 1) (g*)' 



L°°([0,iJ]) 



L°°{[0,T]xI 



we get: 



p'{t)-a\t) 



< te 



Kot 



L°°{[0,i?]) 

TY{pl) + Ca 



L°°([0,i]xl 











1' 


f 1 


Vdivy*(r,x) 


dx dr 









+ • 

L°°([0,/J]) Jo Jr<* 

Using hypothesis (I), we obtain 
p'{t)-a\t) 



div h/^(r,x) -Ty^(r,x) 



dxdr . 



< tCi\\r-s\ 



Li 

L°°([0,t];Li(K'*;R")) 



te 



+e 



Cd 



L°°([0,R]) 

TV(p;) + 



L°°([0,i?.]) 



Vdiv tr*(x) 



+ C/||r| 



L°°([0,i],Li(K<*;R")) 



L°°([0,-R]) 









L°°([0,i?.])_ 



Hence, for T small enough, we can apply the Banach Fixed Point Theorem. 

Now, using the total variation estimate (|4.4p . a standard iteration procedure allows to 
obtain global in time existence. 

Total Variation Estimate. To prove the estimate on the total variation, apply Lemma [4.1l 
and use a procedure entirely similar to that exploited in the proof of Theorem 12.21 
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Stability. We use [TH Proposition 2.10] to obtain 

\\plit) - P2{t)\\^^ 

< ||Po,l -Po,2||li +e'^°*(TV(po,i)+iQlki|lLoo(C7+ llVdivlTillLi)) 



^ ^1 



L°°{[0,R]) 



f X(pi(r))-X(p2(r)) 
Jo 



dr 



+t\\q[\ 



L-{[0,«])l 



""2 - + + II ^^2|lL°°)|kl - 121 



L°°([0,iJ]) 



dr 



+ lki|lL-([o,i?]) divX(pi(r)) -divX(p2(-r)) 
+*lki|lL-([o,i?])l|div ^2 -div ITiIIli + t(C/ + ||div t^2|lLi)lki - 92|lL-([o,i?]) 

< ||Po,l — Po,2||li 

(<^/ + ll^2|lL-)e''°* 



X 



TV(po,i) +tCrf||gi||L^([o^^])(C7 + llVdivlTillLi)) \\q'i - 92||Loo([o,ij]) 

+ (C/ + ||div t^2|lLi)lkl - 92|lL-{[0,iJ]) 

e"°1k'i||Loo([o,/j]) (TV(/5o,i)+tCrf||gi||L^(fo_^])(C/ + ||Vdiv7i|lLi 



) ) II ^^2 - t^l 



+ lkl|lL-([0,R])l|div(tr2 - iTl)! 



(TV(po,i)+tQ||<?i||Loo([o,^j)(C,+ ||Vdiv7l||Li) 



ll''l|lL°°([0,fl]) 

+ \\Qi\\l^{[o,r]) 

X /" IIpiC"^) -P2(t)||li dr . 
Jo 

Let us denote 

t + ||72|lL-)e''°* (tV(po,i) + tCrf||gi||Loo([o,fi]) (C7/ + || V div 7i ||li 



a{t) 



X\\(ll ^2||L«={[0,i?]) 
+ lldiv l^2|lLi)lkl -92|lL-([0,i?]) + lklllL-([0,RJ) 

^TV + t Qllgi ||loo([o,«]) (C, + II V div 7i IIli 



||div (1^2 - 



Ip'^ot II II 

||yi|iL°°([o,fi]) 



x||'y2 - t^l|lL°° 



lL°°{[0,i?]) 



b{t) = Ci^ 

+ ll'?lllL°°([0,iJ]) 

so that, by integration, we get 

||piW-P2(t)||Li < (l + te*''^ 



TV(/5o,i) + tCrf||gi||Loo([o,R])(C/ + ||VdivTri||Li) 



|Po,l - /5o,2||ti + , 
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which completes the proof. 



□ 
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